Growing networks under geographical constraints 
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Inspired by the structure of technological networks, we discuss network evolution mechanisms 
which give rise to topological properties found in real spatial networks. Thus, the peculiar structure 
of transport and distribution networks is fundamentally determined by two factors. These are the 
dependence of the spatial interaction range of vertices on the vertex attractiveness (or importance 
within the network) and on the inhomogeneous distribution of vertices in space. We propose 
and analyse numerically a simple model based on these generating mechanisms which seems, for 
instance, to be able to reproduce known structural features of the Internet. 

PACS numbers: 05.50.+q, 89.75.Hc 



43 

Of 

6 
o 

CO 
CO 

O 

CO 

>v 

43 ■ 

Oh: 



> 

o 
o 

o 
c- 
o 

CO 

o 

CO ' 
>>'. 

43 

Or 



X 



Technological networks, such as transportation or com- 
munication networks, are man-made networks designed 
for transport of resources between sites distributed over 
a certain geographical area [J . Depending on the type of 
network, the resources can be information, wares, elec- 
tricity, or persons, and the geographical area range from 
a small region to the whole world. Examples of techno- 
logical networks are, among others, the Internet 0, ||, 
the airline @, i, H and railway @, S] networks, and the 
electric power grid [J, Q . The Internet consists of a set 
of routers linked by optical fibre or other type of physi- 
cal connection, and it turned into an indispensable tool 
to get information from and about whatever part of the 
world. The airline network, which principal function is to 
transport persons and wares, has all the airports of the 
world as vertices of the network, and the corresponding 
nonstop scheduled flights connecting the airports as its 
edges. Electric power grids, on the other hand, are sets 
of generators, transformers, or substations connected by 
high-voltage transmission lines. 

The most prominent feature of these technological net- 
works is that they are embbeded in a real physical space, 
with vertices having well-defined positions. This is not 
the case of other types of networks, such as citation or 
biochemical networks, in which the pos ition of vertices 
has no physical meaning (see [H, OS lU El for general 
reviews). In many communication and transportation 
networks, the cost of establishing long-range connections 
between distant spots is usually higher than the cost of 
establishing short-range connections. This is clear for 
networks such as the Internet or the railway network, 
where establishing a long-range conection is obviously 
expensive because long channels need a larger infrastruc- 
ture. For electric power grids, the connection cost be- 
tween farther spots is even higher, given that in long 
high-voltage lines a large amount of energy is lost during 
the transmission. 

This dependence of the connection cost on the distance 
is one of the most prominent mechanisms governing the 



evolution of technological networks, and it is determinant 
for understanding their structure. As a consequence of 
that, for example, not all connections between vertices 
are equally probable; neighbouring vertices tend to con- 
nect to each other with larger probability than distant 
ones. This, in turn, is the origin of some of their more 
characteristic properties. 

The most important quantities designed for capturing 
network's structure are the degree distribution, the av- 
erage distance between the vertices, and the mean and 
local clustering coefficients. The degree distribution P(fc) 
gives the probability that a randomly selected vertex of 
the network has degree k, i. e., that it is connected to 
k other different vertices. Most technological networks 
exhibit degree distributions that decay as a power-law 
P(k) ~ fc -7 , i.e., they exhibit a scale-free character. 
However, power grids or railway networks, in which long- 
range connections practically do not exist, typically show 
exponential degree distributions [f| Q . The average path 
length I is defined as the mean distance between each 
two vertices in the network, where the distance between 
any two vertices is defined as the number of edges along 
the shortest path connecting them. Finally, the cluster- 
ing coefficients measure the local tendency of vertices to 
form highly connected clusters. The clustering coefficient 
of a vertex is defined as the ratio between the number of 
connections existing among its nearest topological neigh- 
bours -the vertices which are connected through an edge 
with it- and the maximal number of edges which can exist 
among them. The mean and the local clustering coeffi- 
cients, C and C(k), are the averages of the clustering 
coefficients over all vertices of the network, and over all 
vertices of degree fc, respectively. 

Large mean clustering coefficients and average path 
lengths of many technological networks can be under- 
stood taking into account their growth mechanisms. 
Thus, the fact that vertices tend to link to their "phys- 
ical" neighbours yields a large probability that, given a 
vertex in the network, its "topological" neighbours are 
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also connected between themselves. That typically gives 
rise to large values of mean clustering coefficient. On 
the other hand, since distant vertices tend to be poorly 
connected between themselves, shortest paths connecting 
farther nodes are usually long, and pass through many 
vertices in between. Statistical measures of the lengths of 
connections confirm that the large part of edges in most 
transportation and distribution networks are short-range 
connections @. 

The distance between the vertices is an important pa- 
rameter in the modelling of transportation networks, but 
the "cost-distance" dependence is not the only mech- 
anism responsible for their structure. In general, the 
structure of technological networks is both a function 
of what is geographically feasible and what is techno- 
logically desirable. For reasons of efficiency, some long- 
range connections are typically always present, in spite 
of their high cost: in many cases, connecting two dis- 
tant vertices trough a long necklace of neighbouring ver- 
tices slows down the global transport in the network and 
makes it inefficient. Long-range connections are observed 
both in the Internet and the Airline network [f|. Ad- 
ditionally, when long-range connections exist, they usu- 
ally link the highly connected vertices of the network 
0, E|. That is not surprising. If a telecommunication 
company or an airline decides to make a big investment 
in creating a long-range transport channel, it typically 
wants to link sites which are somehow important and 
well-connected (depending on the type of network, they 
can be technological, touristic or commercial spots), so 
that the amount of information or wares which will be 
exchanged between them compensates the expense. 

Networks embedded in a metric space with distance- 
dependent connection probabilities are called spatial or 
geographical networks [l2| . In the past few years several 
models have been proposed in order to study their struc- 
tura]prmertie^JIpJl3J^ 
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Most of them combine the preferential attachment mech- 
anism [lH , which is widely accepted as the probable ex- 
planation for power law degree distributions seen in many 
networks, and distance effects. The last typically lead to 
a deviation from the scale-free behaviour when the dis- 
tance constraints are sufficiently strong. Almost all these 
studies have focused on the effects of geography on the 
degree distribution, ignoring other important character- 
istics. These are, however, of primary interest. Thus, 
in a exhaustive study of the Internet [|[ Vazquez and 
collaborators found not only that it is a scale-free net- 
work, but also that the local clustering coefficient C(k) 
and the nearest neighbours' average function k nn fall as 
power-law functions with exponents —0.75 and —0.5, re- 
spectively. On the other hand, Gastncr and Newman 
|6| showed that strong geographical constraints tend to 
produce networks with an effective network dimension d 
close to d ~ 2. These new quantities are essential when 
measuring such important features as degree-degree cor- 
relations and the hierarchy Q and planar [|| characters 



of networks. 

In this paper we propose several network evolving 
mechanisms which are able to develop spatial networks 
that exhibit most of the features found in real technologi- 
cal systems, i.e. reproducing correct values for C(k), k nn 
and d. The basic properties of our generating principles 
are: 

i) The knowledge of any given vertex about the net- 
work is limited to a certain (Euclidean) neighbourhood 
of the vertex (a property of locality). Each vertex is 
"aware" of the characteristics of all vertices belonging to 
its neighbourhood, but not of the characteristics of the 
rest of the vertices of the network. 

ii) The range of this physical neighbourhood is gov- 
erned by a cost function which establishes the importance 
of the geographical constraints. As the connection cost 
grows, the range of the neighbourhood decreases. 

iii) As usual, the network grows by adding vertices and 
edges. At each time step, new vertices are added and 
connected to the system; additionally, new edges may be 
set between vertices already existing in the network. 

iv) Preferential attachment condition. Vertices try to 
connect to vertices of large degree - the more attractive 
ones - lying within their neighbourhood. 

Apart from these requirements (already considered e.g. 
in Ref. [l||), we add two new ingredients: 

v) The interaction cost governing the range of each 
neighbourhood depends on the attractiveness of the ver- 
tex associated; the larger is the vertex attractiveness the 
larger its interaction range. 

vi) The probability that a new vertex appears in an 
isolated area, geographically far from the rest of the ver- 
tices, is smaller than the probability that the new vertex 
appears close to an already existing vertex. (Regarding 
to this last point, a similar idea has recently been pro- 
posed by Kaiser and Hilgetag [H[). In addition, vertices 
cannot appear too close to each other. 

The last two conditions are inspired by the properties 
of technological networks. While condition v) is obvi- 
ously a consequence of the fact that the degree of im- 
portant vertices in the network grow faster the "richer" 
they already are (preferential attachment), there is also 
another reason. It is also because they extend their "ten- 
tacles" more far away. Condition vi) mirrors the fact that 
vertices do not appear over a geographical area at ran- 
dom. Consider, for instance, the Internet and the elec- 
tric power grid. When a power plant is constructed in a 
region too far from civilisation, the plant supplies elec- 
tricity to the buildings close to it, but no high- voltage 
transmission lines link the station to the grid of the 
civilised world if the distance is large; the plant usually 
remains isolated. (In fact, electric power plants are not 
constructed far from civilisation. The inhabitants of iso- 
lated regions usually use small generators for personal 
use. The station is constructed only when civilisation 
comes to the region.) In the Internet, routers concen- 
trate in towns, rather than in deserted areas. That is 
quite natural, people live and work in towns! Conse- 
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quently, new Internet accesses tend to appear in towns, 
in the vicinity of other already existing accesses. Fur- 
thermore, the more industrialised is the town, the more 
rapidly the number of Internet accesses grows. On the 
other hand, we must also take into account, that vertices 
do not appear extremely close to each other. Thus, con- 
structing two big power stations in close vicinity to one 
another (for example, one kilometre apart) is not rea- 



sonable; it is cheaper to build one bigger station which 
supplies electricity to the entire area. It also is not com- 
mon for a family house to have two routers, since one 
router can supply Internet to all computers in the house. 
Therefore, we assume that there are certain areas, not too 
far but also not too close to the existing vertices, where 
new vertices will more probably arise than in others. 




FIG. 1: (ai), (&i), and (ci): graphical representations of the models (a), (b), and (c), respectively. The three networks have 
1000 vertices; their size is L = 1881 (ai), L = 1982 (02), and L = 1982 (03). Note that all edges of network (ai) are short-range 
connections, while in network (ci) edges connecting distant vertices do exist. (02), (fo), and (C2): Degree distribution of the 
networks represented in (ai), (02), and (03), respectively. Models (b) and (c) are scale-free. The slope of the straight lines are 
-2.24 (62) and -2.95 (C2). (03), (fe), and (03): Local clustering coefficient C(k) of the networks represented in (ai), (^2), and 
(03), respectively. The behaviour of C(k) follows power laws for three models. The slope of the straight lines are -0.40 (03), 
-0.46 (63) and -0.58 (&3). Notice the double logarithmic scales in all graphs. 
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These basic ideas can be implemented in very different 
ways giving rise to different growth models. Consider, for 
instance, the preferential attachment prescription. One 
has to decide whether the attachment probability de- 
pends linearl y o n the degree, as in the Barabasi- Albert 
construction [38[ , or whether it must follow another law. 
One also has to decide about the interaction's depen- 
dence on distance and on vertex attractiveness. More- 
over, the probability of the appearance of a new vertex 
at position r can be a complex function involving the po- 
sitions of all already existing vertices in the network, or 
depending only on the position of the vertices closer to 
the point. 

The appropriate implementation of a determined pat- 
tern depends obviously on the particular geographical 
system that we attempt to model. Numerical simulations 
show that different (but similarly oriented) prescriptions 
produce models showing qualitatively the same behav- 
ior. This fact supports the general value of the principles 
proposed. Since we are not yet intending to studying any 
particular real-world network, but are interested only in 
capturing some general features of spatial systems, we 
will adopt here the simplest realization of the guiding 
prescriptions. 

We start from a preselected area in a two-dimensional 
Euclidean plane. In this area, we place at random m 
vertices, so that the distance between any two of these 
initial vertices is larger than a given r m i n , the minimum 
distance that will separate vertices in the network. Since 
the m vertices are placed at random, the order of mag- 
nitud of separation between them depends on the size 
of the preselected area. Now we let our network grow 
around these initial vertices. At each time step a new 
vertex with mi proper links is added to the network and 
connected to mi vertices already present. Additionally, 
once the new vertex is attached, new edges are dis- 
tributed among all the vertices of the network. In both 
cases, vertices and edges are added to the network only 
if the geographical constraints allow for the addition. 

Our geographic constraints are defined by the two char- 
acteristic distances r TO j„ and r max , which define a ring 
area around a point. At each step we choose at random 
a vertex of the preexisting network. From this point, us- 
ing polar coordinates, we put a new vertex at a position 
given by a radius r and an angle <j) picked up at random 
from homogeneous distributions r m i n < r < r max and 
<f> G ]0, 211]. If this new vertex happens to be at a dis- 
tance smaller than r m i n from some preexisting one, the 
selection is rejected and a different old vertex is chosen. 
Note that this prescription does not give an homogeneous 
distribution in space when r max — > oo and r m i n — > 0, but 
cssencially means that smaller distancies to the chooscn 
vertex are preferred. 

In order to connect the new vertex to the system, we 
consider the nodes of the network within the circle of ra- 
dius r max from the newly introduced one n. This circular 
area around the new node is considered to be its physical 
neighborhood. If the number of vertices in the neighbor- 



hood is smaller than mi, the newly introduced vertex is 
connected to all of them; if their numer exceeds mi than 
it is connected to exactly the mi ones with higher degree. 
Note that the fact that the range of the neighbourhood 
is precisely r max ensures that a new node is connected to 
at least one old one. 

The second process, consisting of the addition of new 
edges between vertices, works in a similar way. We ran- 
domly choose a vertex v of the network, and then, from 
the vertices that belong to its physical neighbourhood 
but are not yet connected to it, we choose the mi ver- 
tices having larger degree and connect v to them. Here, 
however, the interaction range r v of the neighbourhood 
of v is governed by the function 

Tv ^max ~t~ Pk v , ( 1 ) 

where k v is the degree of vertex v, and (3 and 7 are non- 
negative tuning parameters whose function is basically 
to define the area which a vertex "sees" depending on its 
importance (degree) within the network. In case that the 
number of vertices that belong to the neighbourhood of 
v, and that are not yet connected to it, is q < TO2, then 
only q edges are added to the network. Note that the 
effects of geography disappear when r max — > 00. 

Let us comment on some aspects of the model. First, 
vertices are not distributed at random over the area of 
study, but their distribution depends on the "history" 
of the network. New vertices appear in the vicinity of 
the vertices already present in the network. Second, the 
interaction range of a vertex is a function of its attractive- 
ness, or importance in the network. If a vertex increases 
its importance in the network, then its interaction range 
grows. (Here, we do not take into account the fact that 
old vertices can remain obsolete with time.) Third, we 
impose the following preferential attachment condition: 
connect to the more attractive vertex of the network that 
you can "see" . 

Extensive numerical simulations confirm that this sim- 
ple model is able to reproduce many of the properties of 
spatial networks. In the present study we restrict our- 
selves to two sets of values for the parameters of the 
model. The first one, which includes three different spa- 
tial cases, illustrates the impact of the cost-distance di- 
chotomy on network structure. We consider the follow- 
ing values: mi = 1, rri2 = 1, r m i n = 500 m.u., and 
Tmax = 1000 m.u. (where m.u. stands for an arbitrary 
"metric unit"). The fact that we impose r max = 2 ■ r m i n , 
i. e., that r max is only twice r m i n , indicates that in this 
case we deal with networks for which the cost of new 
vertices establishing long-range connections is very high. 
Additionally, we choose mo = 7, and a radius of 14000 
m.u. for our initially preselected disc-area, within which 
the mo = 7 initial vertices may be placed at random. 
The three cases we distinguish are: Case (a), j3 = 1 and 
7 = 1.4, corresponding to a spatial network in which the 
geographical constraints are extremely important (in this 
case long-range connections are practically inexistent); 
Case (b), j3 = 1.5 and 7 = 2.3, an intermediate case; 
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FIG. 2: 77(C) as a function of the L. Fr om bottom to top, 
models a, b, and c. The black straight line correspond to a 
network of dimension 2. The results are used to estimate the 
dimension of the three spatial networks considered (see text 
for more details) . Notice the double logarithmic scales of the 
picture. 



Case (c), j3 = 2 and 7 = 4, for which vertices of high 
degree are allowed to establish long-range connections. 
The selected values of parameters are certainly arbitrary 
and are adopted in order to illustrate the effects of the 
distance-cost-dependence. In effect, for (3 = 0, the result- 
ing network is practically a tree, since no edges can be 
placed between old vertices, while for very large f3 and 7 
a "winner-takes-all" phenomenon emerges, in which al- 
most all vertices are connected to one super-hub with an 
enormous degree. 

Figure [T] compares the results of simulations corre- 
sponding to these three cases. To be able to draw the 
resulting networks, we consider small graphs with only 
1000 vertices (numerical simulations indicate that the 
structure does not significantly change as the order of 
the networks grows). Panels a\, bi, and c\ show the ef- 
fects of the selective growth of the interaction range with 
the degree of vertices: For systems where long-range con- 
nections are highly expensive (model a), even the most 
important vertices of the network are connected only to 
a few close neighbours. As the cost of establishing long- 
range interactions decreases, connections between distant 
vertices in the network begin to appear, in particular, be- 
tween high degree vertices (models b and c). 

The degree distribution evidently changes as the ge- 
ographical constraints are gradually loosened. Thus, 
model a shows a degree distribution which decays ap- 
proximately exponentially (panel 02); no vertices of high 
degree can be found. The degree distributions of mod- 
els b and c, however, exhibit well-defined power law tails 
(in spite of the small order of the networks considered): 
P{k) ~ fc- 2 - 25 (panel b 2 ) and P{k) ~ fc~ 2 - 95 (panel c 2 ). 
The corresponding local clustering coefficients also show 



FIG. 3: Layer average degree, (k)c, as a function of shell 
number C. The curve with the highest peak corresponds to 
model c, the intermediate one and the flat one to models b, 
and a, respectively. 



a behaviour very close to that found in real networks 
(see panels 03, 63, and C3). All three models exhibit 
power law behaviours for C(k): C(k) ~ fc~ - 40 (panel 03), 
C(k) ~ fc" 46 (panel 63), and C(k) ~ fc" 58 (panel c 3 ). 
In addition, the mean clustering coefficient C of these 
spatial models is always quite large, about C ~ 0.33 
for all three of them. The number of triangles (cycles 
of length three) in the network is, however, 668 (model 
a), 1593 (model b), and 1341 (model c). On the other 
hand, the average path length decreases as the amount 
of long-range connections grows from I = 20.18 (model 
a)to I = 8.83 (model b) and I = 5.01 (model c). This 
result is quite natural, and shows the transition from a 
quasi-planar graph with a structure quite similar to a 
lattice (model a) to a typical complex network structure 
found in most geographical networks (models b and c). 

Making use of the effective dimension d, Gastner 
and Newman showed that, networks where geograph- 
ical effects are extreme are essentially planar graphs 
(i.e., they can be drawn on a map without any edges 
crossing). The effective dimension can be defined as 
d = lim£^ 00 logA7'(£)/log£, where A/"(£) is the average 
number of vertices which can be found within a distance 
of C steps or less from a vertex. In finite networks no 
limit C — > 00 can be taken, but good results for d can 
be achieved by plotting log M(C) against log£ for the 
central vertices of the network and measuring the slope 
of the resulting line (of course, far away of the saturation 
region corresponding to exhausting of the network). Cen- 
tral vertices are those vertices of the network that have 
minimum eccentricity, being defined as the maximum dis- 
tance from the vertex to any other vertex in the network. 
(Note that central vertices are sometimes defined as the 
vertices having larger "betweenness centrality" , as in 0] . 
We use, however, the classical definition from the graph 
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FIG. 4: Analysis of the degree-degree correlations of the 
three models considered. Nearest neighbours' average func- 
tion k n n{j) against j. From bottom to top, models a, b, and 
c. Note that models b and c exhibit dissortative mixing, while 
model a is slightly assortative. 



theory.) 

Figure [2] shows on double logarithmic scales how 77(C) 
behaves as a function of C. From bottom to top the 
curves correspond to model (a), the critical behaviour 
AT(jC) ~ C 2 , model (b), and model (c). We see that the 
effective dimension of model a is certainly smaller than 
two, which is not a surprize provided that model a creates 
practically a planar graph. The dimensions of models b 
and c -which are obviously not planar graphs- are larger 
than two. The difference between the three models at 
Af(C = 1) is due to the fact that central vertices are 
usually the more connected ones of the network. 

Tomographic studies reveal interesting details too. To- 
mography deals with the study of the structure of lay- 
ers which surround a given vertex (the root) in the 
network [3^, 0, Hl|. The principal motivation for ex- 
amining the tomography of a network results from its 
importance for understanding the spreading phenomena 
taking place in networks. We concentrate here on the 
layer average degree (k)c = J2kkPc(k), where Pc, the 
degree distribution in shell C, is defined as Pc(k) = 
(Er N cA k ))/(E k ,r N cAk)), with NcAk) being the 
number of vertices of degree k in layer C for root r. The 
study of (k)c for the three networks considered shows a 
peak whose height decreases as the cost of establishing 
long-range connections grows (fig. [3]). The results indi- 
cate that the mean degree (k)c=i increases rapidly as the 
number of long-range connections and hubs in the net- 
work grows. On the other hand, the average shell degree 
decreases rapidly for more distant layers, C > 1. This 
interesting result shows that vertices with large degrees 
are rapidly exhausted in this type of networks, which has 
especial importance when dealing with spreading phe- 
nomena, like spreading of information or infections. Note 



that this result has important effects on epidemiological 
properties: vertices of large degree are rapidly affected 
by the spreading of an infection. On the other hand, 
in a network like that of model (a) the propagation of 
any spreading agent will be similar to the propagation 
on a lattice: the spreading agent will primarily reach the 
nearest physical neighbours. 

In figure [1] we show the correlation properties of these 
three models. Degree-degree correlations are determined 
by the probability function which gives the proba- 
bility that a randomly selected edge connects one ver- 
tex of degree % to another of degree j. Thus, a net- 
work is said to be degree-degree uncorrelated if £y = 
(2 — Sij)iP(i)jP(j)/(i) 2 , which only means that the prob- 
ability that an edge connects to a vertex of a certain de- 
gree k is independent from whatever vertex is attached 
to the other end if the edge. Otherwise, the network is 
said to be degree-degree correlated. Most real networks 
are correlated, and usually exhibit either "assortative" 
or "dissortative" mixing [42|. Assortativity means that 
high-degree vertices attach preferably to other highly 
connected vertices, i. e., with a larger probability than 
in uncorrelated networks; on the other hand, dissorta- 
tivity stands for when high-degree vertices tend to con- 
nect to low-degree vertices, and vice versa. Thus, a very 
useful quantity for measuring the correlation's degree of 
a network is the nearest neighbours' average function 
k nn (j), which expressed in terms of £y, can be writ- 
ten as ~k~nn(j) = *(! + <5i.j)£tf)/(X)i(l + hj)£ij). It 

takes the constant value k nn (j) = (j 2 )/(j) if no type of 
degree-degree correlation exist, while it is a decreasing 
(increasing) function if dissortative (assortative) mixing 
is present. In the picture we plot k nn (J) as a function of 
j. The lowest curve, corresponding to model (a), shows 
then that the network is slightly assortative. This feature 
of model (a) is due to the fact that the areas containing 
a large density of vertices usually contain a large density 
of edges (see Figure [TJ ai), corresponding probably to 
important areas of the space; and vice versa, the areas 
containing a small density of vertices also contain a small 
density of edges. On the other hand, models (b) and (c), 
in which the geographical constraints are not so strong, 
present dissortative mixing. Interestingly, for both mod- 
els k nn (j) falls with j following power laws of the form 
knn(j) ~ j~ L , just like happens in real networks. 

Model (a) reproduces quite well the properties of those 
systems where vertices and edges are embedded in the 
two-dimensional physical space, like, for example, electric 
power grids or road networks. However, none of the three 
models considered above is suitable for characterising 
world-scale systems such as the Internet or the network of 
airline routes. The reason is that, in such large-scale sys- 
tems, vertices are usually not uniformly distributed in the 
region under study (as occurs in our preceding models, 
see pictures [T]ai; b\, and ci), but they concentrate in a 
number of technological areas distributed over the world. 
Thus, a more realistic model for describing such systems 
must take into account that in large-scale geographical 
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FIG. 5: (a): Graphical representation of a small network (N = 1000 and L = 1985) corresponding to model (d). Note that 
vertices concentrate in certain areas of the space, and that the long-range connections of the network link end- vertices of large 
degree, (b): Degree distribution of model (d). (c): Degree-dependent clustering coefficient C(k) of the model, (d): nearest 
neighbours' average function k nn as function of the degree j. Note that P(k), C(k), and k nn fall off as power law functions, 
(e): Average degree, (k)c, as a function of shell number C. From (k)c=i the average degree decays interestingly following a 
power law (inset of the picture). 



networks there usually exist many "desert" regions lying 
between the areas where vertices can be found in abun- 
dace. Such a pattern is easy to construct by varying the 
ratio between r max and r min in our model. This aspect 



is actually considered by our second selection of parame- 
ters. As we will see next, inhomogeneous distribution of 
vertices in space influences quantitatively the statistical 
properties of networks. 
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Let us thus consider a larger ratio for r max /r m i n , for 
example, r max = 5 • r TO j„. (Note that the change in 
the ratio r max /r m i n modifies not only the distribution 
of vertices in space, but also makes cheaper the cost of 
establishing connections for new vertices, since the neigh- 
bourhood of a new vertex will, in comparison, contain 
more vertices.) The values of the parameters for this 
model (model d) are now the following: too = 7, m\ = 1, 
m 2 = 1, r min = 200 m. u., r max = 1000 m. u., j3 = 2 
and 7 = 3. As before, the initial area, where the too = 7 
vertices are randomly placed, has a radius of 14000 m. 
u. and the order of the network is N = 1000. Figure [5] 
(a) shows this network in the two-dimensional Euclidean 
space: The model simulates perfectly the tendency of ver- 
tices to concentrate in different areas having a high den- 
sity of vertices (as if these areas were urban centres, i.e., 
cities, or city agglomerations), which are linked through 
long-range connections which join vertices of large de- 
gree, usually belonging to different geographical commu- 
nities. 

The properties of this last construction are especially 
interesting, since they reproduce many structural proper- 
ties found in the Internet. Thus, the degree distribution 
of the model follows a power law P(k) ~ k~ 2A2 (panel 
b of figure [5]). The mean clustering coefficient is large. 
For our small network C ~ 0.7. The dcgrcc-dcpcndcnt 
clustering coefficient C(k) decays as a power law function 
C(k) ~ fc" - 97 (panel c of the Figure). The decay of C(k) 
gets however slower for networks of larger size N. The av- 
erage path length is very small, I = 3.74 for the parame- 
ters as in Figure O and numerical simulations with larger 
networks indicate a small-world behaviour. In addition, 
the network shows dissortative mixing: the nearest neigh- 
bours' average function k nn decreases as k nn {j) ~ j~ - 52 



(see Fig. \5\d). The coincidence of these properties with 
the ones of the Internet is astonishing: Let us now remind 
that i) the degree distribution of the Internet follows a 
power law with exponent 7 ~ —2.5, ii) the local cluster- 
ing function C(k) behaves as C(k) ~ k~ ' 75 , and iii) k nn 
decreases with j following the function k nn {j) ~ j -0 ' 5 [![• 
Finally, in Figure[S]we plot (e) the average degree {k) c as 
a function of shell number C, corresponding to the study 
of tomography. We see again that hubs are found only a 
few steps away from any vertex, and interestingly, that 
(k)c drops as a perfect power law from C = 1 on (see 
inset of the Figure; note the double logarithmic scales). 

To conclude, we introduce several network-generating 
mechanisms taking into account the constraints that ge- 
ography impose on the evolution of large-scale network 
systems in physical space. We suggest that two proper- 
ties are determinant for the structure of such geographi- 
cal networks: the fact that the spatial interaction range of 
vertices depends on the vertex attractiveness and the fact 
that that vertices are not randomly distributed in space. 
Simple implementations of these mechanisms show that 
the essential difference between "strong geographical" 
networks, such as electric power grids, and "weak geo- 
graphical" networks, such as the Internet or the airline 
network, could be the cost (economical or technologi- 
cal) of establishing long-range connections. On the other 
hand, inhomogeneous distribution of vertices in large- 
scale networks seems certainly to be a relevant generat- 
ing element of their hierarchical character. In any case, 
the agreement of our results with the properties found 
in real networks suggest that the mechanisms proposed 
may play a key role in the evolution and structure of 
networks. 
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